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ABSTRACT 

In this report, the basic heat transfer and isothermal flow friction 
characteristics for five perforated nickel matrices, one perforated nickel 
fin with solid splitter matrix, and one solid nickel matrix, all of 
"parallel plate" configuration are presented. In addition, similar in* 
formation is provided for a 20° skew matrix fabricated from perforated 
nickel. Finally, performance of two offset rectangular fin matrices 
constructed from stainless steel, is evaluated. Where meaningful, a 
comparison of the performance of these surfaces is made. All matrices 
are investigated by the transient method of testing utilizing the maximum 
slope analysis of Locke and include the effects of longitudinal thermal 
conduction as reported by Howard. 
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NOMENCLATURE 



English Letter Symbols 

2 

A Matrix total heat transfer area, ft 

2 

A £ Matrix minimum free flow area, ft 

A^ r Matrix total frontal area, ft^ 

A g Solid matrix cross sectional area available for thermal 

conduction, ft^ 

2 

A^ Conduction area corrected for effect of perforations, ft 

A* Plane surface area, ft^ 

a Plate thickness, ft 

a Short side of a rectangular flow passage, ft 

b Plate spacing, ft or in. 

b Long side of a rectangular flow passage, ft 

C Flow stream capacity rate, (mc^), Btu (hr 0 F) 

C Matrix capacity, W c , Btu/ °F 

S S 

c Specific heat at constant pressure, Btu/(lb °F) 

P 

c s Matrix material specific heat, Btu/(lb °F) 

Djj Hydraulic diameter, 4r^ 

d Inside pipe diameter, in. 

» 

2 

E Flow friction power per unit area, HP/ft 

f Mean friction factor, dimensionless. This is the "small" 

or "Fanning" friction factor. (Ratio of wall shear stress 
to the fluid dynamic head.) 

G Exchanger flow stream mass velocity, (m/A ), lb/(hr ft^) 

c 

g Proportionality factor in Newton's Second Law, g = 32.2 

c (lb/#) (ft/sec^) 

h Unit conductance for thermal convection heat transfer, 

Btu/ (hr ft^ °F), or heat transfer power per unit area 
per degree temperature difference, Btu/(hr ft^ °F) 
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English Letter Symbols (continued) 



j 


2/3 

Colburn factor = N Np^ , heat transfer characteristic, 

dimensionless 


J' 


Colburn j factor based on plane surface area, A*, dimension- 
less 


K c ,K e 


Contraction loss coefficient for flow at heat exchanger 
entrance or exit respectively, dimensionless 


k 


Unit thermal conductivity, Btu/(hr ft^ °F/ft) 


k s 


Matrix thermal conductivity, Btu/(hr ft °F/ft) 


L 


Total matrix flow length, ft 


• 

m 


Mass flow rate, lb/hr 


p 


2 

Pressure, #/ft 


p 


Porosity for matrix surfaces, dimensionless 


q 


Heat transfer rate, Btu/hr 


R 


Universal gas constant, ft-#/lb- 0 R, (53.35 for air) 


r h 


Hydraulic radius, (A C L/A), ft, (4r^ = hydraulic diameter) 


s 


Solidity of a perforated material, 1-3>J dimensionless 


T 


Absolute temperature, degrees Rankine, °R 


t 


Temperature, degrees Fahrenheit, °F 


V m 


Matrix volume, ft** 


V 

s 


Material volume corrected for effects of perforations 


W 

s 


Mass of matrix, lb 


X 


Distance along the flow passage in direction of flow, ft 



Greek Letter Symbols 



oi 


Aspect ratio of a rectangular flow passage, b/a, dimension- 
less 


<3 


Compactness; ratio of total heat transfer area to the volume 
between the plates, ft^/ft 3 


A 


Denotes difference 


8 


Fin t'hir.knpfi.q. ft 



Fin thickness, ft 



x 



Greek Letter Symbols (continued) 



* 


Conduction area reduction ratio due to perforations in 
perforated material; cross sectional area of perforations/ 
cross sectional area, dimensionless 


e 


Time 


\ 


Longitudinal conduction parameter, k A , dimensionless 

s s 

mLc 

P 


r 


Time parameter 




Fluid viscosity coefficient, lb/hr-ft 


y 

e 


Area reduction ratio due to perforations in perforated 
material, dimensionless 

Density, lb/ft^ 


cr 


Ratio of freeflow area to frontal area, A c /Af r , dimension- 
less 


P 


Denotes "function of" 


5 

I 


Compactness for perforated material including the effect 

of area reduction, A*s/V , dimensionless 
J nr 

Summation 



Subscripts 



atm 


Local atmosphere 


f 


Fluid 


i 


Initial, individual 


k 


Equivalent 


m 


Mean or matrix, as appropriate 


o 


At orifice 


s 


(Matrix material) solid 


STD 


Standard temperature and pressure 


1 


Upstream 


2 


Downstream 



xi 



Force and Mass Units 



lb 


Denotes pounds mass to differentiate from 


# 


Denoting pounds force 


Dimensionless 


Groupings 


n r 


Reynolds number, (4r^G/yu), a flow modulus 


n R(p) 


Reynolds number for pipe, (dG/yu) 


N st 


Stanton number, (h/Gc ), a heat transfer modulus 

P 


%u 


Nusselt number, (h4r^/k), a heat transfer modulus 


N Pr 


Prandtl number, (yiiCp/k), a fluid properties modulus 


N Tu 


Number of heat transfer units, (hA/mCp) 


j 


Generalized heat transfer grouping - Colburn "j n factor, 
(N 8t N p This factor versus N^ defines heat transfer 

characteristics of the surface. 


f 


Mean friction factor. This is the "small" or M Fanning n 
friction factor (Ratio of wall shear stress to the fluid 
dynamic head). This factor versus N^ defines the fric- 
tion characteristics of the surface. 




Conduction parameter, k A /mLcp for solid material; 

k A. /mLc for perforated material 
s k p 




Equivalent conduction parameter (corrected for equivalent 
length in perforated material) L 


r 


Time parameter, hA g 

W c 
s s 
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INTRODUCTION 



While it is appreciated that the designer of a compact heat exchanger 
is confronted with many qualitative evaluations, nevertheless he is general 
ly working specifically within limitations relating to flow frontal area, 
weight, pressure drop and volume. It is for this reason that different 
heat transfer surfaces may be effectively compared on graphs using co- 
ordinates of f and j (or f/j ratio) versus N„. 

As in many other disciplines, model theory is the most satisfactory 
approach from a cost effectiveness standpoint* The use of such methods 
as steady state "steam to air" tests or "transient" testing techniques 
provides the basic flow friction and heat transfer characteristics without 
necessitating the construction of a full scale heat exchanger* 

The data presented herein was obtained by the use of a transient test 
technique or "single blow". C7] This technique is one of heating a sample 
core to some predetermined value and then after thermal equilibrium is 
achieved in the specimen, instantaneously removing the upstream heat source 
and monitoring the downstream temperature response of the fluid at the core 
exit. 

The objectives of this thesis specifically were: to modify an experi- 

mental test rig to operate with a larger capacity blower and to test sever- 
al offset rectangular fin, of stainless steel, and "modified parallel" plate 
of perforated nickel, sample cores. 

In the testing of perforated nickel specimens, the effect of varying 
shape and size of perforations was to be investigated. 
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SUMMARY OF THEORY 



The "starting stage of operation", "single blow", or "transient test 
technique", as it has been called by various authors, had its inception 
with the analytical treatment by Anzelius in 1926. Nusselt in 1927, Hausen 
in 1927 and 1929, and finally Schumann in 1929 treated the problem of ana- 
lytically determining the fluid and solid temperature as a function of 
position and time for a porous solid, initially at a uniform temperature, 
when subjected to a change in the entering fluid temperature. Some repeti- 
tion occurred as the authors were unaware of the work performed by their 
contemporaries at the time C9] « Anzelius and Schumann made similar as- 
sumptions; however, the results are different. Anzelius simplified the 
problem; whereas, Schumann proceeded to obtain an explicit mathematical 
solution for the temperature distribution. This solution was for a liquid 
flowing through a porous solid and consisted of two infinite series which 
were related to Bessel functions. Schumann states Cl7j that if the fluid 
is a gas, the solution to the problem is more formidable. He proceeds to 
relate how dimensional analysis may be used so that his investigation re- 
mains usable for compressible fluids. 

The porous solid is considered to be homogeneous and at a uniform 
temperature throughout with the fluid flowing at this same temperature. 

With a sudden change in the entering fluid temperature, the temperature 
of the fluid and the solid as a function of both position and time re- 
mains to be evaluated. The theoretical generalized heating curves are 
developed from an energy balance on the porous solid as presented in [19j 
and restated in Cll . Inasmuch as experimental procedure used in this 
report is dependent upon the assumptions and boundary conditions of this 
energy balance, it is briefly summarized below: 
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Assumptions made in the analysis are: 

a. Properties of the fluid are temperature independent 

b. Fluid flow is steady 

c. Porous solid is homogeneous 

d. Thermal conductivity of both fluid and solid is infinite 
perpendicular to the fluid flow direction 

e. Thermal conductivity of fluid is zero in the flow direction 
Boundary conditions are Cl3} : 

(a) The matrix is initially at a uniform temperature 

(b) At time equal zero, the temperature of the entering fluid 
changes instantaneously to a different, constant value, 
i.e., a step change in fluid temperature 

(c) The matrix boundaries are adiabatic 

A I 

Heat absorbed by solid = ^s^$ c s dx 

Heat transferred to the solid by convection * hbC^f _ ^-s^ y 

Heat transferred from the fluid by convection = mCe^P d* 

dx 

Heat transferredfrom solid by conduction = -k s A s dx 

d* 2 - 

Thus an energy balance for the fluid and solid yields: 

(1) (^ s A s c s dx = k s As^sdx + hb 

d© 

(2) mcp |£f dx + hb (-t f -t s )dlx =0 
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Rearranging the above equations and substituting in the following para- 



meters [3] 



T' - generalized time variable 



= LA 

W s c 



s'-s 



©- 

m 



t] 



where h = unit conductance for convective heat transfer, 

Btu/hr ft 2 °F 

2 

A = matrix heat transfer area, ft 
W g c g =* matrix capacity, Btu/ °F 
© = time, hrs. 

= mass of fluid contained within the matrix, lb. 
m = mass flow rate of fluid, lb/hr. 

x = distance from entrance of matrix along flow passage 
in the direction of flow, ft. 

L = matrix flow length, ft. 

Regrouping the above equation for / T > and introducing the specific heat of 
the fluid, c p , we get: ^ « _hA_ q _ UW f C p * 

W S C S m W s c s c ? L 

Since the thermal capacity of the fluid contained at any time within the 

matrix (W f c ) is much less than the thermal capacity of the matrix (W c ), 
I p s s 

the second term can be neglected; accordingly: 



LA 

W s c s 



0 



(3) 



z = generalized position variable = 



Inasmuch as the number of transfer units (N_ ) = LA- 

1 Cp 

becomes 



h A x 
•^Cp L 

, this simply 



z = 




( 4 ) 



4 



( 5 ) 



\ ■ k A 
v _8 8 

me L 
P 



2 

where k^ ■ matrix thermal conductivity, Btu/hr ft °F / ft, 

A = solid matrix cross sectional area available for thermal conduc- 
8 tion, ft^. 



m, c , and L are as given above. 



Equations (1) and (2), together with (3), (4), and (5) yield: 



dr 



X N 



TU 







( 6 ) 






where t = temperature of the solid 
s 

tf = temperature of the fluid 

For the case of thermal conduction in the solid parallel to the direction 
of flow equal to zero, equations (6) and (7) simplify: 







( 8 ) 



= -fc -tr (9) 

T s r 

dz 



It is with these two differential equations, that Schumann utilizing 

Bessel functions determined the time - temperature history for the solid 

after a step change in fluid temperature. For a more thorough treatise, 

the reader is referred to references [13] or Cl7] . 

The task of determining the N value for a test point required the 

Tu 

superposition of the experimentally obtained time - temperature heating 
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or cooling curve on the family of theoretical curves. Although this would 
uniquely determine the value, it can be appreciated as being a tedious 
effort. 

Fortunately, in 1950, Locke [13] developed a "maximum slope" technique 
which bypasses this task. He differentiated the theoretical curves and 
plotted the slope of the curve versus 'T' /N^ u with N as a parameter. The 
expression for the slope is 

d [VN Tu ] ^ Kxu - r 



-T,(Zvfn^r)e' (HT ^ Y) 



where the downstream fluid temperature is evaluated at x = L so that z = 

N_ and t, = t-„. 

Tu f f2 

In the efforts of both Locke and Schumann, the thermal conductivity 
of the solid parallel to the flow was taken as zero. Although the effect 
of longitudinal conduction is not pronounced at high Reynolds number range, 
it is considerable at low laminar Reynolds numbers. Several authors have 
considered the effects of longitudinal conduction, introducing the conduc- 
tion parameter, A. WhenAl^ O , the maximum slope is not a unique 
function of N^ u but depends on X as well. Howard [7] has tabulated maxi- 
mum slopes as functions of and A. , as well as presenting graphs, 
which he evaluated by a finite difference technique with the aid of a 
digital computer. The evaluation of in this report is based on these 
tabulated values in so far as practicable. 

Satisfying the aforementioned idealizations, initial and boundary 
conditions, leaves the determination of the conduction parameter, A- , 
and the experimental evaluation of the maximum slope of the heating or 
cooling curve as the only two remaining items necessary to determine the 
heat transfer characteristics of the matrix. 
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EXPERIMENTAL TECHNIQUES 



The results of Howard's finite difference calculations in reference 

[7] were used in this report for evaluating N and ultimately the 

Tu 

Colburn j factor. For this reason, the idealizations imposed by Howard 
were restrictions placed on the experimental set up used. Howard £7] 
considered transient heating of a porous solid when subjected to a step 
change in fluid temperature including the material longitudinal conduc- 
tion effects. The idealizations are restated below: 

"(1) The flow entering the matrix is steady and uniform in 

velocity and temperature, and remains steady and uniform 
at any cross-section as it flows through the matrix. The 
thermal conductivity of the matrix is infinite in the 
direction normal to the flow and finite in the direction 
parallel to the flow. Thus the problem is one-dimensional 
in space. 

(2) The thermal capacity of the fluid contained at anytime 
within the matrix is small compared with the thermal 
capacity of the matrix. That is, the fluid is normally 
restricted to a gas and the results to be presented 

. will not be expected to apply for liquid fluids . This 
also means that for the fluid there will be no time de- 
pendent terms in the equations. 

(3) The thermal properties of the fluid and matrix are con- 
stant and uniform. 

(4) The convective heat transfer coefficient is some suitable 
average and remains uniform and constant. 

(5) At time zero the change in fluid temperature will be a 
step function, with the matrix and its entrained fluid 
initially at some uniform and equal temperature." 

Idealization (1) was achieved, as far as steadiness and uniformity of 
temperature and velocity were concerned, by responsible experimentation. 

In an effort to improve on the previously used test apparatus of Bannon 
[1] , a new inlet cone was designed and the test section relocated. The 
inlet cone was fashioned according to the specifications of Smith and 
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Wang [l8] for contracting cones giving uniform velocity at the throat. 
Although this design was specifically for circular cones, the resulting 
square inlet cone fabricated, together with wire screen type flow straight- 
eners provided the required uniform velocity profile. Velocity profiles, 
taken at the matrix inlet over the entire flow range, verified required 
uniformity. The test section was placed upstream of the air metering 
device so that full advantage of the uniform velocity profile could be 
realized. Uniform temperature profiles within +0.5 °F were obtained by a 
wire screen type heater explained in detail in APPENDIX B. The latter part 
of the first idealization pertaining to the thermal conductivity of the 
matrix cannot be achieved precisely; however, the effects of deviating 
from it can be made negligible in the solid, if resistance to heat transfer 
is small within the solid compared to the heat transfer resistance be- 
tween solid and gas. In other words, "transverse conduction" effects have 
been neglected; whereas, the effects of longitudinal conduction are consid- 
ered with the use of X. > the conduction parameter. In a matrix where 
adequate mixing takes place, no criterion is necessary for "transverse 
conduction" effects. 

The second idealization is self explanatory, i.e., the fluid is re- 
stricted to a ga$. 

Idealization (3) for the thermal properties can be guaranteed if the 
temperature differences are kept small. A maximum air temperature varia- 
tion from 70 ° F to 90 0 F results in a maximum deviation from the arith- 
metic average in the order of + 1.57* in viscosity and is negligible in 
specific heat. 

The convective heat transfer coefficient presented as the final result 
of this experiment is an overall quantity; accordingly, idealization (4) 
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is satisfied 



Satisfying the fifth idealization of a step change in fluid tempera- 
ture is perhaps the most rigorous requirement of all. Several investiga- 
tors have attempted to meet this requirement by various schemes. The 
technique of Bannon £l} and Theoclitus and Eckrich {j20j appears to be 
most satisfactory; therefore, it was the one selected for this investiga- 
tion. 

Figure 1 displays a schematic of the experimental set up. Figure 2, 
is a photograph of the complete test apparatus. The system is comprised 
of an air source, flow straightener screen section, heater section, matrix 
test section, air metering device, and pressure and temperature measuring 
systems. Each of these are discussed at length in APPENDIX B. 

Air was drawn through the "wind tunnel" by a Spencer Turboblower. 

The 0.0031" diameter electric wire heaters upstream of the matrix were 
energized so as to heat the air to roughly 20 degrees above ambient. The 
upstream and downstream fluid temperatures at the matrix test section were 
monitored to ascertain when the matrix had been heated to a uniform tempera- 
ture. When this condition was satisfied, the heaters were quickly shut 
down and the time-temperature history of the entering and exit air at the 
matrix was recorded on a Minneapolis Honeywell "Brown" Recorder. The 
recorded curve is the experimental cooling curve from which the maximum 
slope was determined. 

The mass rate of flow of air was measured by an ASME orifice meter 
with pressure taps one diameter upstream and one - half diameter down- 
stream, using a range of orifice plates. 

Flow friction performance of the matrix was obtained by measuring the 
isothermal pressure drop across the matrix and the mass rate of flow of 
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air through the matrix. Pressure taps upstream of the matrix and upstream 
of the orifice meter were used to record static pressure entering the 
system and at the orifice section respectively. 



Data reduction relations are discussed in APPENDIX A. 
calculation for one test run is presented in APPENDIX D. 
The following data was recorded for each test run: 



(a) 


Orifice diameter (d ) 

o 






(b) 


Entering static pressure 


(P s> 


(c) 


Pressure drop across 


the 


matrix ( A P m ) 


(d) 


Static pressure upstream 


of orifice (P ) 
o 


(e) 


Pressure drop across 


the 


orifice (AP ) 
o 


(f) 


Orifice temperature 


< fc o> 




(g) 


Atmospheric pressure 


( p „ ) 
atm 


(h) 


Time - temperature cooling curve 


(i) 


Recorder chart speed 







A sample 
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4. DESCRIPTION OF SURFACES 

Surfaces examined in this experiment fell into one of two major classi- 
fications - perforated nickel plate matrices or offset rectangular fin of 
stainless steel* In the former case, two of the specimens examined - 
the 20 ° skew matrix and a modified parallel plate matrix, both of 160/40TV 
perforated nickel - had previously been examined by Bannon [l] . These parti- 
cular surfaces were selected to provide an additional check on the system 
reliability* 

The 20 ° skew matrix was formed by corrugated sheets, which were skewed 
relative to one another to prevent meshing of adjacent sheets at an angle 
of 20° . This matrix configuration unlike both the modified parallel 
plate and the offset rectangular fin matrices does not require splitter 
plates* 

The "modified parallel plate" configuration consists of formed sheets 
separated by splitter plates to prevent nesting. The assembled matrix is 
a maze of rectangular passages with a value of oC , the aspect ratio of a 
rectangular flow passage, = 7. The following sketch depicts the ideally 
formed "modified parallel plate" matrix: 
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A modified parallel plate matrix of each type of perforated nickel was 
constructed and tested to examine the effects of varying shape, quantity and 
size of perforations on this configuration. 

The offset rectangular fin surfaces were provided by the courtesy of 
the AiRsearch Manufacturing Company, Los Angeles. The method of identifi- 
cation employed with this configuration is best illustrated by the follow- 
ing example: 2 j| 

16R - .153 - .143 (0) - .004 - SS 

<2)0 (I) 0 © © © 

Number of fins per inch 

(2) Type of fin flow cross-section (R = rectangular) 

( 5 ) Fin height, inches 

@ Uninterrupted fin length in the flow direction, inches 

® Type of surface ( 0 = offset) 

( 6 ) Fin metal thickness, £ , inches 

@ Fin material (SS = stainless steel, A1 - Aluminum) 

Two matrices of this configuration were investigated: - the 16R - 
.153 - .143 (0) - .004 - SS and the 20R - . 100 - .125 (0) - .004 - SS. 

The offset - rectangular fin matrix is comprised of offset fins, separated 
by splitter plates (0.004" thick) to prevent nesting of the fins which 
would alter the flow characteristics. The isometric sketch below best 
identifies this plate fin arrangement: 
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5. 



PRESENTATION OF RESULTS 



Heat transfer and flow friction data for each matrix investigated are 
exhibited in both tabular and graphical form, TABLES I through X and 
Figures 9 through 18, using the Colburn j modulus, Fanning friction factor, 
and Reynolds number. In evaluating the Reynolds number, the hydraulic 
diameter, which is defined as four times the hydraulic radius, was used. 

The effects of entrance, exit, and flow acceleration have been considered 
in the evaluation of f. (See Appendix A). 

Geometric and physical properties of the surfaces are given in Figures 
3 through 8. Figures 19 through 23 show heat transfer and friction compari- 
sons between various matrices investigated. In Figures 24 and 25, M figure 
of merit" type curves are presented. The j/f versus Reynolds number thus 
plotted gives an indication of the required matrix flow frontal area for 
a given pressure drop. 

Lastly, in Figures 26 and 27, heat transfer power versus flow friction 
power curves on a unit area basis, evaluated for fluid properties at stand- 
ard conditions of dry air at 500 °F and one atmosphere, are represented. 

(See Appendix A). 
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6 . 



EXPERIMENTAL UNCERTAINTIES 



Those items of possible error arising from lack of exact duplication 
of the idealizations, boundary and initial conditions are discussed in 
the "Experimental Technique” portion of this report and no further attempt 
was made to evaluate this type of error. 

The experimental errors which do not fall into the above area are 
three fold and are discussed below. These three categories are: (1) 

uncertainty of physical constants, (2) inaccuracy in the geometric meas- 
urements and (3) instrumentation inaccuracies. For the determination of 
the uncertainties reported, the method of Kline and McClintock jjL2j was 
used, 

(1) Values for the physical constants required were obtained 
from [4j , and [6^j . The uncertainties in these values, as best 

as can be determined, are listed below: 



c 


+ 


0.5% 


s 






c 


+ 


0.5% 


p 






N Pr 


+ 


2.0% 


k 


+ 


0.5% 


s 






JU 


+ 


1.0% 



It is noteworthy that the accuracy for the value of c g depends upon 
the metal used. Where pure metals are used the accuracy quoted is feasible; 
however, this is not the case with alloys. Thus the tolerance of c g can be 
a major cause for inaccuracy ^13^ • 

(2) Due to inconsistencies in construction of the matrices, 
errors in lineal dimensions are considered to be + 0.5%, Inasmuch as the 
weight of the matrix can be determined as accurately as desired, the error 
in W g is considered negligible. Based on this, the uncertainties in geo- 
metric measurements are considered to be as follows: 
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+ 1.07. 



A, A fr> A^, A g> A k> 

L + 0.5% 

W < 0.1%, negligible 

s 

(3) Instrumentation inaccuracies were essentially those 
obtained in pressure measurement. With the exception of the temperature 
measurement at the orifice, all temperatures recorded were temperature 
differences which were recorded in inches (See Appendix A). The tempera- 
ture measurement at the orifice was obtained by a copper-constantan shield 
ed thermocouple which was read by a Leeds and Northrup portable potentio- 
meter. Assuming adequate manufacturer’s calibration of the thermocouple 
wire, the estimated possible error is + one-half of the smallest division 
on the potentiometer, or + 0.0025 millivolts which is approximately +0.1 
°F. 

Since the range of pressures varied, requiring several different mano 
meters, the maximum uncertainty encountered was used in the uncertainty 
analysis. These values are listed below: 



p„ 


+ 


1.25% 


o 




A*, 


+ 


, 1.25% 


P 


+ 


1.70% 


m 






P _ 

atm 


+ 


0.0005" Hg (negligible) 



One final parameter for which uncertainty must be determined before 

an analysis of the experimental results can be made is . Inasmuch as 

Colburn j factor is not a linear function of Reynolds number, uncertainty 

for high and low N values was considered. Using Figure 30 and the un- 
it 

certainty in Maximum Slope as 2.0% these values were determinable. 

^ . dN Tu 

For N_ =3.0, A- = 0 m =+7.2% 

Tu w Tu 
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For N = 25.0, A. = 0.05 

Tu * 



dN Tu 



N. 



= + 10.07. 



Tu 



Now, with the aid of reference M , the uncertainties in j, m, N^, and f 
were evaluated based on 20: 1 odds. 






( wA 






_Vbv l / 


v jw,. / 


v ^v n ; j 



■»& 



where w , w„, ....w = uncertainty interval based on certain odds, 

1 * 2 n J 

and w^ is the interval for the result. 

R 

Considering the evaluation of the Colburn j factor: 



then: 



2/3 



j = N N 

Tu A Pr 



Wj = 






'/? 



JfT. 

(|lir^AcN P r /3 



Dividing by j yields: 



= 

j 



^Nm \ 
v Nru / 



^ ( Wk,\\ (-L vJ Np ^ 
A / + l Ac / W Kpr " 1 



i K. 
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For the N_ 
Tu 



X = 0 case: 



= 3.0, 



Wj _ 



z x 

(.002-)% (.01) Hoi) C^) J " 



For the N = 25.0, X = 0.05 case: 

Tu 



^ = j\.i) +(.0i)\O)\ ('f) 



= .1018 



a 7. 5 % 

j _ 



*» = \O.Z% 

j 



By similar analysis the uncertainties at 20:1 odds were determined 
for the following: 



m + 1.07. 

N r + 2.37. 

f + 4.37. 
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7. 



DISCUSSION OF RESULTS 



In order to verify the reliability of the data obtained by the test 
apparatus used, several surfaces on which data was currently available 
were tested. Two of these matrices, the 160/40TV perforated nickel 
M parallel plate 11 and 20° skew, were evaluated by Bannon £l] . These two 
matrices were originally evaluated by the same transient test technique 
on basically the same test apparatus; therefore, satisfactory comparison 
of the data provided the necessary confidence in the test rig. Figures 
9 and 16, as well as TABLES I and VIII, present the results. 

The two offset rectangular fin matrices tested were essentially tested 
by Briggs in [ 2 ] with the exception of slightly thicker splitter material 
and using aluminum in lieu of the stainless steel. The data in [2^ was ob- 
tained by a steady state "steam- to-air n test technique. Although this 
technique permitted investigation atf much higher values of Reynolds number, 
a sufficient overlap in the Reynold’s range existed so that a comparison 
was possible. A good comparison for the 20R -.100 - .125 - .004 - SS 
surface for Colburn j factor exists; however, the friction factor re- 
ported in [ 2 ] differs noticeably. Results are shown in Figure 17 and 
TABLE IX. For the 16R - .153 - .143 - .004 -SS surface, the values for 
Colburn j factor for Reynolds number < 1000 are in excellent agreement 
with those presented in £ 2 ] . In the transition range, the j values re- 
ported herein are noticeably different; however, the values for friction 
factor agree nicely throughout the entire duplicated Reynolds number range. 
It is important to note here that by the "maximum slope" method of analysis 
an inherently large error is possible in the low N region (1.5 ^ N Tu “ 
3.5) M . The size of this error in N relative to an error in maximum 
slope can best be appreciated from Figure 30. Consider, for example, a 
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27. error in maximum slope with an N value of 2.3 and A_ ~ 0. The 

Tu 

percent error in N , from Figure 30, is of the order of 207.. For values 
of between 2.0 and 2.3, the percent error has not been evaluated. The 
1.5 — N^ u ^ 3.5 region coincides with the Reynolds range for which the 
reported j values differ from [2], and may account for the difference 
noted. The results for this matrix appear in Figure 18 and TABLE X. 

In spite of the areas of disagreement, the correlation was such that 
the data for the j in the Reynolds range 1000 and the friction factor 
data for these matrices are considered reliable. It is noteworthy that 
the surface tested in [ 2 ] is specified as 15.76R - .153 - .143 -.004 - Al, 
so that the matrix evaluated herein may only be an approximate duplicate. 

A comparison of the relative performance for these matrices is presented 
in Figure 23. Additionally, Figure 25 illustrates j/f vs. N R , a flow 
area "goodness" factor. 

The perforated nickel matrices investigated covered a wide range of 
percent open area, 12 to 50%. Both slotted and round perforated nickel 
surfaces were evaluated for the "modified parallel plate" geometry. The 
perforated nickel geometric properties are given in APPENDIX C for the 
various types investigated. In order to make a reasonable comparison of 
the performance among these matrices, the same frontal area and flow length 
were specified. Results are presented in Figures 9 and 10 and listed in 
TABLES I and II for the slotted perforated nickel matrices; those for the 
round perforated matrices are given in Figures 11, 12, and 13 and tabu- 
lated in TABLES III, IV, and V. 

In order to establish a meaningful fiducial point, a solid nickel 
matrix of the identical material (obtained from the same manufacturer) 
with the same frontal area and flow length as the perforated matrices was 
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fabricated. The performance of this matrix is illustrated in Figure 15 
with tabular values given in TABLE VII. 

The best performing slotted perforated nickel matrix was the 160/40Q. 

The 50G round perforated nickel (50% open area) was the best surface in- 
vestigated among the ’’parallel plate" nickel matrices. 

As another approach to evaluating the Colburn j factor, j* was con- 
sidered. The parameter j 1 is based on the plane surface area, A*, as if 
no perforations were present; whereas, j is based on the heat transfer sur- 
face area, A, which includes the solidity correction, s. The value of j' 
is readily attainable from the relationship 

j 1 = j x s , where s = A/ A*. (See Appendix D) . 

The determination of the solidity is fully described in APPENDIX C. Figure 
13 for the 50G perforated nickel matrix includes the values for j 1 versus 
Reynold's number in addition to the f and j values. 

A comparison of the performance of the two slotted perforated matrices 
together with the solid nickel matrix is shown in Figure 19. A similar 
comparison for the round perforated nickel matrices is given in Figure 20. 

It is readily apparent that the perforated material yields a higher 
heat transfer characteristic. Inasmuch as the friction factor also shows 
an increase, a more appropriate comparison can be made from the flow area 
"goodness" factor curve in Figure 24. The heat transfer power (h ) versus 

olll 

flow friction power (E CTn ) per unit area given in Figures 26 and 27, for the 
slotted and round perforated nickel matrices respectively, is another "good- 
ness" factor evaluation. Here the information for each of the matrices is 
presented based on standard conditions for fluid properties and a common 
hydraulic diameter. (See APPENDIX A). 
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So that an appreciation for the perforated material might be realized 
in direct transfer type of heat exchanger application in addition to its 
use in matrix type for rotary regenerators, a "parallel plate" matrix with 
the 160/40TV perforated nickel fins and solid nickel splitter plates was 
fabricated. The performance of this matrix is given in graphical and 
tabular form in Figure 14 and TABLE VI. A comparison with the solid 
nickel matrix is made in Figure 21. Figure 24 contains the j/f presenta- 
tion and Figure 26 shows the ti _ vs. E for this matrix. Here again is 

a noticeable increase in heat transfer performance, but with a negligible 
increase in friction factor. 

For the 160/40TV perforated nickel material, a 20° skew matrix was 
compared in performance against the "modified parallel plate" to see how 
these two configurations compare with one another. The separate perform- 
ance for each has been mentioned previously; however, the direct comparison 
is shown in Figure 22. As with the other nickel surfaces, the flow area 
"goodness" factor is represented in Figure 24. 

In the results of this investigation are several items which are general 
in nature, but nonetheless extremely important. In the laminar Reynolds 
number range, friction factor is inversely proportional to Reynolds number. 
By Reynolds analogy, it was anticipated that the Colburn j factor should 
likewise be inversely proportional to Reynolds number in this region. 

However, the experimentally obtained values of j in the low laminar Reynolds 
range utilizing the equivalent conduction length (defined in APPENDIX C) 
together with TABLE XI or Figure 29 yielded lower values for j. For this 
reason, a re-evaluation of the conduction parameter, , and maximum slope 
data points for Figure 29 was performed. The equivalent conduction para- 
meter X was computed and the value was determined by utilization of 

the extrapolated j vs. N curve, since j is proportional to N . This new 

JK 1U 
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determined value of N was considered to be the actual N for the parti- 

Tu T u 

cular "maximum slope" and conduction parameter. Values determined in this 
fashion are asterisked in TABLES I through VI and VIII. The tabulated 
conduction parameter values in these tables are the equivalent conduction 
parameter (See APPENDIX D). 

The second item is that of an equivalent conduction length which was 
quoted earlier. In the perforated fin material the slots or round holes, 
as the case may be, provide for a winding heat flow path. It was appreci- 
ated that the matrix flow length, L, was not the true conduction path but 
somewhat less than actual. A mean conduction path length evaluation is 
presented in APPENDIX C. The conduction parameter, /V , based on the flow 

length, L, was evaluated, and by simply multiplying by the ratio L/L (also 

k 

given in APPENDIX C), the value of X was determined. 

k 

A third item is the inherent error in N in the range 1.5^ N m — 3.5 

Tu Tu 

which was mentioned earlier. Howard \j7 J published an error curve for N , 

Tu 

including slope vs. N , see Figure 30. 

Tu 

Fourthly, the entrance and exit length effects for a matrix with a 
small hydraulic radius are small. In the evaluation for Fanning friction 
factor (APPENDIX A), the effects of exit and entrance length were consider- 
ed using values of K and K from Figures 37, 38, or 39, as appropriate. 

For the perforated surfaces, the K values were taken corresponding to 
laminar flow conditions; however, in the offset rectangular fin geometry 

the values were taken from those Figures for N = 00 in accordance with 

R 

reference [ii] • The offset fins never have velocity profiles that are 
fully established, that being a purpose in the design. The same calcula- 
tions were repeated on the offset rectangular fins using the K values for 
laminar flow. The friction curves plotted for the offset rectangular fins 
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were coincident with those obtained earlier; however, some change was noted 
in tabulated values. 

A final item worth mentioning is again concerned with the low Reynolds 

range. A slight alteration in the conduction parameter value has a marked 

effect on the value in this range, which is apparent from Figure 29. 

For this reason, when utilizing a nickel material, or any other material 

with a wide range of specified values for thermal conductivity, it is 

mandatory that the appropriate value of k be obtained for calculating 

s 

purposes. The value of A. is directly proportional to this k g value, 
i.e. , 

/V = k s A s 
L m c p 
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CONCLUSIONS 



Reliable data is obtainable from the experimental apparatus based on 
a comparison with those surfaces previously evaluated elsewhere. 

A very definite improvement in performance is achieved with the utiliza- 
tion of a perforated material. Kays, in reference £l0] > postulated the 
improved heat transfer characteristic without a "large amount of form drag 
so characteristic of high performance surfaces". The perforations apparent- 
ly disturb the heat transfer sub boundary layer so that tubulence increases 
and the resulting high convection conductance is realized. Inasmuch as the 
perforations are fine, the downstream edge is considered to be in the wake 
of the separated region. The perforated nickel matrices investigated here- 
in confirm Kays* hypothesis. Unfortunately, inasmuch as the perforated sur- 
face investigated in £l0} is greatly different in geometric and physical 
properties, no meaningful comparison with those tested in this investiga- 
tion could be made. 

As stated previously, the experimentally obtained Colburn j values in 

the low laminar Reynolds range using TABLE XI or Figure 29 were lower than 

j values obtained by extrapolating the j versus N curve. Inasmuch as fric- 

R 

tion factor and Colburn j factor are both inversely proportional to Reynolds 
number in the laminar region owing to Reynolds analogy, the values for 
Colburn j factor extrapolated from the experimental curves are considered 
more reliable in the low laminar range. 

Effects of longitudinal conduction have been considered in the analysis 
of data. For the perforated material matrices, corrections for the actual 
conduction path length and for the area removed by the perforations have 
been included. 

The best performance for a perforated "parallel plate" matrix was at- 
tained by the 50G type material. For the offset rectangular fin surfaces, 
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the 20R - .100 - .125 - .004 - SS matrix was superior. 
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9. 



RECOMMENDATIONS FOR FURTHER STUDY 



1. Due to the inherent error in the high Reynolds number range where 
1.5 ^ N f^3.5, a more reliable means for determining N^^ values should 
be investigated. Perhaps pursuit of Hausen's cyclic operation , a 
steady state periodic-flow technique, is appropriate in that this method 
measures effectiveness from which N^ can be determined subsequently. 

Since the results obtained using Schumann's analysis, even with longitudin- 
al conduction effects considered, at the low Reynold's range appears in- 
adequate; the use of another testing method may prove fruitful. 

2. Inasmuch as the best performing surface for perforated nickel 
was the 50G type (50% open area), other types, such as the 40G which has 
a slightly larger hole diameter, with same plate thickness and percent 
open area as the 50G may prove worthy of test. 
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Figure 1. Schematic Diagram of Test Apparatus 
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Figure 4« Geometric and Physical Properties of Round Perforated 

Nickel ( 125 M, 125 P, and 50 G ) "Parallel Plate" Matrices 
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Matrix Material 


Solid Nickel 


Specific Heat (c s ) Btu/lb °F 


0.1065 


Thermal Conductivity (* s ) 3tu/ hr ft °F 


38.7 


Material Thickness, incncs 


0.0020 


Total Heat Transfer Area (A) f t^ 


20.1875 


Frontal Area (Af r ) ft/ 


0.0695i 


Total Conduction *rea (A s ) ft/ 


0.010 u 


Free Flow area (A c ) ft 


O.o5939 


a 

Matrix Volume. (V m ) ft v 


0.011588 


Matrix Density ( m ) ib/ft^ 


77.7 


Hydraulic Diameter (D^) ft 


0.001961 


Compactness (0) ft^/ft^ 


■1742 .1 


Porosity (p) 


0.854 



Figure 6. Geometric and Physical Properties of a 
Solid hickel ”P fir all ei Plate” Matrix 






Matrix Material 


20 Degree Skew 
160/40 TV 
Perf .Nickel 


Specific Heat (c s ) Btu/lb F 


0.1065 


Thermal Conductivity (k s ) Btu/ hr ft F 


38.7 


Material Thickness, inches 


0.0022 


Total Heat Transfer Area (a) ft^ 


15.41 


Frontal Area (Af r ) ft^ 


0.0734 


Total Conduction Area (A^) ft^ 


0.00349 


Free Flow Area (A c ) ft^ 


0.06412 


Matrix Volume (V m ) ft^ 


0.01211 


Matrix Density (^ m ) lb/ft ^ 


50.1 


Hydraulic Diameter (D^) ft 


0.00252 


Compactness ( Q ) ft 2/ft- 


1272.8 


Porosity (p) 


0.873 



Figure 7. Geometric and Physical Properties of a Slotted 
Perforated foickel (160/40 TV) 20° Skew Matrix 




I 



> , i J i ) )! 

! I 



till 


MM 


lll| i 


IN! 


Mil 


ll 111 


t 1 


Mil 


nm 

1 j 1 1 | 


i ( i 
' ^ \ \ 

i i « 


1 1 1 
i ^ * > 

( 1 


II ! II 

1 1 M 1 


i I 


M ' * 


• ■ • r c 


'j ) ) } 

1 ! ! i 


MM 

1 I j 


) ) pim i 

i 


f 1 1 

>789 
JETT CO. 0 


1 I 1 hi f 1 IT 

1234^6789 

ATHOL.MASS.U.S.A. 



7 8 9 
ATHOLMASS.USA £ 



Matrix Material 


20ti -.100 - 
.J25-.004-SS 


16R -.153 - 
.143-.C04-SS 


Specific Heat (c s ) Btu/lb F 


0.120 


0.120 


Thermal Conductivity (k s ) Btu/ hr ft F 


9.3 


9.3 


Material Thickness, inches 


0.004 


0.004 


Total Heat Transfer Area (A) ft^ 


11.625 


9.0625 


Frontal Area (A^ r ) ft2 


0.07083 


0.C7083 


Total Conduction Area (A s ) ft^ 


0.C1549 


0.01201 


Free Flow Area (A c ) ft^ 


0.05535 


0.05882 


Matrix Volume (V m ) ft^ 


0.0177 


0.0177 


Matrix Density ( ^ m ) lb/ft^ 


66.3 


50.4 


Hydraulic Diameter (Dpj) ft 


0.00476 


0.00649 


Compactness ( (^) ft^/ft^ 


656.5 


511.8 


Porosity (p) 


0.781 


0.830 



Figure 8. Geometric and Physical Properties of Stainless 

Steel Offset Rectangular Fin (20R-.100- e l25-«004 
-SS and l6R-.153-.H3-.OOt -SS) Matrices 





Figure 9« ' Surface Heat Transfer and Friction Data 

"Parallel Plate", 160/^0 TV Perforated Nickel 
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.Figure 10. Surface Heat Transfer and Friction Data 

"Parallel Plate", 160/40 Q Perforated Nickel 
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Figure 11. Surface Heat Transfer and Friction Data 

"Parallel Plate", 125 M Perforated Nickel 
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Figure 12. Surface Heat Transfer and Friction Data 

"Parallel Plate", 125 P Perforated Nickel 
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Figure 13. Surface Heat Transfer and Friction Data 
"Parallel Plate", 50 G Perforated Nickel 
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Figure 14* Surface Heat Transfer and Friction Data 

"Parallel Plate" , 160/40 TV Fina with Solid 
Nickel Splitter Plates 




Figure 15* Surface Heat Transfer and Friction Data 
"Parallel Plate", 0.002" Solid Nickel 
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Figure 16. Surface Heat Transfer and Friction Data 
20* Skew, 160/40 TV Perforated Nickel 



44 






Figure 17. Surface Heat Transfer and Friction Data 

Offset Rectangular Fin, 20ft -.100 -.125 -.004 - SS 










Figure 18. Surface Heat Transfer and Friction Data 

Offset Rectangular Fin, 16R -.153 -#143 -.004 - SS 
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Figure 19. Comparison of Slotted Perforated vs. Solid 
Nickel "Parallel Plate" Matrices 
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Figure 20. Comparison of Round Perforated vs. Solid 
Nickel "Parallel Plate" Matrices 
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Figure 21* Comparison of Slotted Perforated Fin with Solid 

Nickel Splitters vs. Solid Nickel "Parallel Plate" 
Matrices 
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Figure 22. Comparison of Slotted Perforated Nickel "Parallel 
Plate" vs. Slotted Perforated Nickel 20° Skew 
Matrix 
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Figure 23 • Comparison of Offset Rectangular Fin Matrices 
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Figure Zl\, Flow Area Goodness Factors for the Nickel Surfaces 
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Figure 25. Flow Area Goodness Factors for the Offset 
Rectangular Fin Surfaces 
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Figure 26 « Heat. Transfer Power vs» Flow Friction Power per Unit Area 
for the Slotted Perforated Nickel and Solid Nickel 
■ "Parallel Plate" Matrices 
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Heat Transfer Power vs. Flow Friction Power per Unit Area 
for the Round Perforated Nickel and Solid Nickel "Parallel 
Plate" Matrices 
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FIGURE 28. GENERALIZED COOLING CURVE 






56 



NUMBER OF TRANSFER UNITS VERSUS MAXIMUM SLOPE 
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53 
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Figure 30. Error in N^ u Relative to an Error in M axi mum Slope vs. 





Figure 31. Matrix Holder 
(Upstream View) 











Figure 33. Velocity Profile and Temperature Distribution 

Measurement Device 





Figure 34. Mchrorae .lire Heater 
(Upstre-sra View) 



Figure 35 * Close-Up of Apparatus Test 







Figure 37. Entrance and Exit Pressure Loss Coefficients 
for Triangular Cross Section Matrix with 
Abrupt-Contraction Entrance and Abrupt- 
Expansion Exit 
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Figure 38. Entrance and Exit Pressure Loss Coefficients 
for Parallel Plate Cross Section Matrix with 
Abrupt-Contraction Entrance and Abrupt- 
Expansion Exit 
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Figure 39* Entrance and Exit Pressure Loss Coefficients 

for Multiple - Rectangular Cross Section Matrix 
with Abrupt-Contraction Entrance and AbrUpt- 
Expansion Exit 
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TABLE I SUMMARY OF HEAT TRANSFER AND FRICTION RESULTS: 

160/40 TV PERFORATED NICKEL, MODIFIED PARALLEL PLATE 
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From Extrapolated Values of j versus N^ Curve 
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From Extrapolated Values of j versus No Curve 



TABLE III SUMMARY OF HEAT TRANSFER AND FRICTION RESULTS: 

125 M PERFORATED NICKEL, MODIFIED PARALLEL PLATE 
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From Extrapolated Values of j versus Nr, Curve 
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From Extrapolated Values of j versus N Curve 



TABLE Y SUMMARY OF HEAT TRANSFER AND FRICTION RESULTS: 

50 G PERFORATED NICKEL, MODIFIED PARALLEL PLATE 
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TABLE VI SUMMARY OF HEAT TRANSFER AND FRICTION RESULTS: 

160/40TV PERFORATED NICKEL FIN, SOLID SPLITTERS "PARALLEL PLATE 
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TABLE VII SUMMARY OF HEAT TRANSFER AND FRICTION RESULTS: 
SOLID NICKEL, "PARALLEL PLATE" 
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From Extrapolated Values of j versus Curve 



TABLE IX SUMMARY OF HEAT TRANSFER AND FRICTION RESULTS: 
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TABLE X SUMMARY OF HEAT TRANSFER AND FRICTION RESULTS: 
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TABLE XI Tabulation of N«j> u as a Function of the Maximum Slope of the Generalized 
Cooling Curve at Various Levels of the Material Heat Conduction Parameter 
(/M ) Parallel to the Flow Direction 



APPENDIX A 



DATA REDUCTION RELATIONSHIPS 



Summarized herein are those data reduction relationships which are 
of prime importance in obtaining the heat transfer and flow friction 
characteristics as presented from the raw data taken. Before the data 
on the matrices tested can be reduced, or a comparison of the different 

matrices can be made, it is essential to obtain accurate geometrical dimen- 

l 

sions. The geometrical factors of interest are the porosity, p, the hy- 
draulic diameter, D^, and the area compactness, Q . Determination of 
any of these fixes the third since it may be shown that 



(3 = ±£L 
Dh 



[l4] (A-l) 



The definitions of these terms are: 



(1) 


Free Flow Area u t 

p = Frontal Area = 




(A- 2) 


(2) 


4 x Free Flow Area 

D = Heat Transfer Surface Area = 
H 


4 Ac L 
A 


(A- 3) 


(3) 


Heat Transfer Surface Area * 
V? = Frontal Area x Flow Length 


A 

A^ r L 


(A-4) 



MAXIMUM SLOPE 

As presented earlier, the maximum slope of the generalized cooling 
curve is a unique function of N and X. . 







(A- 5) 
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'T* and N^ u have both been previously defined, but for convenience 
the relationships are restated: 



/V'V ^ A 

I 



W s c s 



© and h A 



nu 



m 



C P 



Therefore: 



'i-/n tu .= (r c ?/vhc % ) e 



and 



^(^/Nru.) ~ l ^ c p/vi €> c s ) d0 =■ (c/c s ) d9 



(A- 6) 



Furthermore: 



"tpj_ ~ "t'l 



+ 1 



[7] 



and 






-% -ti 






(A- 7) 



Combining (A«*6) and (A“7) 
± lz~ ti 



J 

lr/H T u] 



c £ 



c 



1 d [if z - t fl ~| 

de 



max 



(A-8) 



From Figure 28: 



( 1 ) 



l d 

dG 



y 

X 



max 



( 2 > x — = a© » sec 

Chart speed 
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(3) y = d( -fc.p. 2 . “ "t-fl ), inches 

(4) M = ("tfj - "ti, ), inches 
Combining these determined values into (A-8) yields: 



( ~ ~tpz_ t j 

L^fi -ti 



C s 

c 



M 



y # Chart speed (A-9) 
X 



I max 



With this value of maximum slope and = 



XI or Figure 29 to obtain a corresponding N 



Tu 



As 

rnCpL 

value. 



enter TABLE 



MASS RATE OF FLOW : 

The mass rate of flow is calculated from the method given in reference 
jj.6} with necessary rearrangements in order to reduce unnecessary hand 
calculations and to gain the full benefits of Murdock £l5~] : 

m = 359Kd e 2 F a Y ]j AV Q V ' (A-10) 



where: 



K = flow coefficient, including velocity of 

c 

approach 3 with C 3 orifice 

coefficient of discharge (obtained from [l5j ) 



and 



Q - 



3 ratio of orifice diameter to internal 



pipe diameter, d^/d 

d 3 orifice diameter, inches 
o 

F^ 3 thermal expansion factor 
Y 3 expansion factor 

vy P 

3 3 specific weight of fluid flowing = 

2 

with P in #/ft 
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R ■ 53.35 ft - it/lb °R for air 
T ■ t(°F) + 459.7 » degrees Rankine (°R) 



AP q " Pressure drop across the orifice, inches 1^0 

substituting P_ for V , since the magnitude of the local acceleration of 
R T 

8 rav *ty, 8> is taken as being equal to the magnitude of standard accelera- 

Q 

tion of gravity, g £ ; and _ ^ - v for K. 
gives: 



ina - 



359 C d Q 2 Y ^ *-El ' 

1 RT 



(A-ll) 



Inserting the numerical value for R and converting P to units of 
#/in 2 yields: 



rn = 



589.81 C Fa Y 

YT T F r " 




- P C-PSta) 

T 



(A-12) 



From Figure 38, reference (l6l : F a ~1.00 

From Figure . 40b, reference ^16^ : '^#“0.9985^ 1.0 
The pressure, P, in (A-12), is hereafter referred to as P^ and is 
determined by the following relationship: 

P, « (P _ - P / 13.6) 0.4912, #/in 2 (A-13) 

j atm o 

Where P is the atmospheric pressure in inches of mercury, 
atm 

P q is the static pressure upstream of the orifice in 
inches of HjO 

13.6 inches of H 0/ inches of mercury 
2 

2 

0.4912 #/in / inches of mercury 
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The temperature, T, in (A-12) is the air temperature at the orifice 



in degrees Rankine (T ) : 

o 



T 0 = t c + 459.7 



(A- 14) 



Replacing F and Y in (A-12) with their numerical values, the equa- 

Bi 

tion reduces to the working form: 



m = 



589.81 C do Z H APo . T> 3 



1 1 ' (3 4 ‘ 



To 



(A-15) 



REYNOLDS NUMBER 

Reynolds number is defined as: 



■ 



(A-16) 



where G is the mass flow velocity based on the free flow area, A^. 



Gr- 



no 

Ar 



no 



? A -fr 



(A-17) 



Substituting (A-17) into (A-16): 



Nts = 






/* A -Pr 



Th 



(A-18) 



But from (A-l) D /p = 



, therefore 



: » = -aa_LL\ ( A- 



19) 



H ' r " R /*-A f rV@, 

Equation (A-19) shows that for a given matrix and mass flow rate the 
Reynolds number is inversely proportional to the compactness. Substitut- 
ing (A-4) into (A-19) yields the working equation: 

4 L m 



Nr " 



j*- A 
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(A-20) 



FANNING FRICTION FACTOR: 



The following sketch together with equation (A- 21) describes the flow 
system under consideration: £llj 



Flow 



AP - 
Pi 




G 2 v x 
2 gc Pi 



l{, \\\\\S\\\\\ WvS 

Cp;-V \W\\ \X\\\ \\N 

r kwwwwxw wHi 



fvjw V\\ VWVWW^I 



lyvwwwwww'ii 

a b 



2 



4 




entrance flow 

effect acceleration 



o 



A v m 
+ f A c 

core 

friction 



- (1 



2 v 2.~ 
p d - K e ) 177- 

exit 

effect 



(A- 21) 



It is important to note that v & ==* and v fe ps in the above sketch 
as the pressure changes (for gas flow heat exchanger application) from 
section 1 to a and b to 2, respectively, are very small relative to the 
total pressure. Inasmuch as the testing is performed with air at moder- 
ate temperatures and pressures, the perfect gas law is applicable (P - ^RT), 



Substituting v = 



and recalling that — !i- = — £. or 



= A c r M . Ac 

£ -- 4 L A l A 

from equation (A-3), equation (A-21) solved for f for the isothermal 



case becomes: 



f = 




?1+ ? 2 ( _*£ 4. K e ) 

2 



con’t on 
next page 
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(con’t.) 



(.H^K 





(A-22) 



where ^ m = ^ ~ with the upstream values subscripted by 1 

and the downstream subscripted by 2. K c is the entrance coefficient and 
K e the exit coefficient, and both are dependent on porosity, shape of 
the flow cross section and the matrix Reynolds number. Values of and 
K 0 are obtained from Figures 37, 38, or 39, which were reprinted from 
reference £ll^J • 

If we consider an order of magnitude approximation, the first term in 
(A-22) is by far the greatest contributor to friction factor for small 
pressure differentials. 

The approximation, therefore, that -p -p 

*1+ 2- = <» S* \ 

reduces (A-22) to the following 



f = j^gc - (K c + Ke)+-^“ (1 + P 2 -)! r H 

J II* 



(A-23) 



Reference ^14] points out that the relationships given here can be 
recombined in such a manner so as to show that f is proportional to the 
porosity cubed and inversely proportional to the compactness (f c>C p / (3 )• 

COLBURN J FACTOR : 

Colburn j factor is defined as: 

j = N st N-p r ^ N Pr (A- 24) 
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Substituting (A-17) for G and multiplying by A/A yields 

3 - is-A _Ac_ N Pr Vs 



m c 



P A 



(A-25) 



but 



u 



m c, 



therefore: 



= N 



j = 



T U 



N 



Tu 



Np, 



V3 



(A-26) 



Combining equations (A-l) and (A-3): 

Ac _ ? 

A §L 



thus: 



J - 



(A- 2 7) 



It is apparent from equation (A-27) that the Colburn j factor is 
directly proportional to porosity and inversely proportional to compact- 
ness P/^). 



HEAT TRANSFER POWER AND FLOW FRICTION POWER : 

An evaluation of the heat transfer power versus flow friction power 
is of interest in that it is a measure of relative performance. 

The heat transfer power per unit area per degree temperature differ- 
ence is £n] = 



h ' £e 4, J_ j 

/V Pr 2/3 Dm J * 



(A-28) 
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Equation (A-28) with the properties evaluated at standard conditions 
of dry air at 500 F and one atmosphere becomes: 

h STD = °' 02195 (•%; )(N R j ) Bcu/hr fc 2 *F (A-29) 

where c = 0.2477 Btu/lb °F 
P 

/*- = 0.0678 lb/hr ft 

^ = 0.0413 lb/ft 3 

N = 0.671 
Pr 

The flow friction power per unit area is [11] : 



E 



1 

2 gc 




(A-30) 



When equation (A-30) is evaluated at standard conditions ( and ^ 
values given above), the flow friction power per unit area at standard 
conditions is: 

E = 1.11 x 10" 7 (-=r* f HP/ft 2 

STD v Uh 7 \lOOO / 

For comparison purposes the surface geometries were reduced to a 

common hydraulic diameter of D = 2 x 10"^ ft. Figures 26 and 27 show such 

n 

comparisons for the nickel surfaces. 
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APPENDIX B 



DESCRIPTION OF EQUIPMENT 

Heat transfer data was obtained utilizing the "single blow" technique 
which consists briefly of monitoring the fluid temperature response at 
the test matrix exit while subjecting the inlet to the test matrix to a 
step change in fluid temperature. 

Friction factor was determined from pressure drop data obtained from 
static pressure taps located in the test section - one at inlet, the other 
at the exit. 

Necessary equipment to perform such an experiment falls into one of the 
following categories: 

(1) Fluid source 

(2) Flow metering system 

(3) Temperature measuring system 

(4) Pressure measuring system 

(5) Fluid heater system 

(6) Matrix holder and test section casing 

FLUID SOURCE : 

Air, the working fluid, was provided to the test apparatus by connect- 
ing the rig to the inlet of a 30HP, multistage, Spencer Turbo-Compressor 
rated for 550cfm operating on a 220 volt a.c. power supply. 

FLOW METERING SYSTEM : (See Figure 1) 

The flow metering system consisted of an ASME standard orifice section 
constructed for d and d/2 pressure taps in a 3.08" inside diameter metal tube. 
M- A wide flow range was obtained by utilizing thin plate concentric 
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orifices of throat diameter 2.310, 1.971, 1.232, 1.540, 0.462, and 0.308 
inches respectively. 

Control of the pressure drop across the orifice was maintained by a 
gate valve downstream of the apparatus, ahead of the compressor inlet and 
another gate valve on the compressor inlet. The blast gate on the compressor 
discharge was pre-set so that the unit could not be operated beyond the full 
load rating of the turbo. 

This turbo-compressor had many distinct advantages over the one used in 
reference some of the more pertinent being (1) a 550 cfm capacity; 

(2) essentially constant pressure throughout the range of the machine; 

(3) the load decreases as less air is used so that the problem of stalling 
at low flow rates is non-existent. 

TEMPERATURE MEASURING SYSTEM : (See Figure 1) 

The majority of the temperatures were measured with iron-constantan 

thermocouples. These thermocouples use iron for the positive conductor 

and constantan for the negative conductor. Measuring T is a group of five 

thermocouples in series, bound together and insulated from each other by 

teflon tape and inserted in a 1/8" diameter aluminum tube to serve as a 

radiation shield. A small- aperature was cut in the inserted tube facing 

upstream. T^ and T^ are five wire thermocouple grids, connected in series, 

mounted permanently in the test section casing. T is a movable five wire 

3 

series connected thermocouple grid located in the matrix holder so that it 
is adjacent to the downstream side of the matrix. The purpose of connecting 
the thermocouples in series was to magnify the emf output so that the instru- 
ment sensitivity to small temperature changes would be enhanced. With T^ 
bucked against T , it is possible to determine the temperature uniformity 
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across the matrix. versus indicates the difference between upstream 
and downstream temperatures; thus, it is this combination which is used to 
record the rate of change of upstream and downstream temperatures with time, 
resulting from the step change in fluid temperature upstream. This differ- 
ential is recorded on a Minneapolis-Honeywell "Brown" Recorder and serves 
as the primary data for heat transfer evaluation of the matrix. The re- 
corder has variable chart speeds so that acceptable cooling curves for 
maximum slope determination can be generated. This instrument has a span 
adjustment which permits the span to be varied continuosuly from the narrow- 
est to the widest span desired from 0 to 55 millivolts. For the work perform- 
ed herein, the recorder was precalibrated for a 0 - 3 millivolt scale. This 
unit also has adjustable supression and damping adjustment features. The 
damping adjustment provides a filter network for removing a.c. strays, and 
provides correct damping so that neither undershoot nor overshoot in the 
response curves is experienced. 

A copper constantan thermocouple was inserted into the duct just a- 

head of the orifice to measure T , orifice temperature. This temperature 

o 

was read in millivolts on a Rubicon Company Portable Precision Potentio- 
meter, converted to degrees Fahrenheit and recorded for each run. 

A copper constantan thermocouple grid was inserted in place of T^ 
to determine the temperature distribution across the cross section during 
initial check out of the apparatus. This thermocouple arrangement was 
part of a special holder lined with balsa wood to provide a smooth, 
continuous flow passage, which also contained an impact tube used for 
checking the velocity profiles (See Figure 33 ). 



PRESSURE MEASURING SYSTEM : (See Figure 1) 

Pressure taps are located upstream and downstream of the matrix and on 
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either side of the orifice. Each tap is connected to an appropriate mano- 
meter or draft gage via Imperial Company "poly-flo** tubing. Quick closing 
valves were installed at various positions in the line to permit isolation 
of sections and as a safety feature. 

The following instruments were utilized: 

(1) Meriam Instrument Company, 0 - 50 M manometer 

(2) Meriam Instrument Company, -8" to +8" manometer 

(3) Ellison Differential Direct Draft Gages, 0 - 6", 0 - 8". 

(4) Ellison Inclined Draft Gages, 0 - 0.5", 0 - 6". 

(5) E* Vernon Hill and Company Type "C" Micromanometer, 

0 - 1.25”. 

Any one of these instruments or a combination can be used to measure the 
differential pressures of the orifice or matrix or the required static 
pressures - upstream of matrix and upstream of the orifice. Cross checking 
of the various instruments assured reliable operation. 

FLUID HEATER SYSTEM : (See Figure 1) 

The heater section is comprised of 28 nichrome wire heaters of 0.0031** 
diameter. The heater system was designed to elevate the air temperature 
20 ° F above ambient for a mass rate of flow of 1000 lb/hour. The nichrome 
material was selected because it has high resistivity, low thermal conduc- 
tivity and specific heat, and has a very small time constant so that it 
permits one to approach the idealized step function. [i] 

The wire heaters are wound two to a bakelite frame, 1/32** between each 
wire, with 50 and 52 wires respectively to each heater on a frame. These 
heaters are connected in parallel, then via a switch to a variable voltage 
bus. The heater frame and its electrical connections are apparent in 
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Figures 34 and 35. 

A schematic wiring diagram for one of the two heater complexes is shown 
below £l"| ♦ A total of 14 frames are wired in this fashion, thereby put- 
ting all 28 heaters in parallel, permitting two to be switched on or off 
at a time. 




The number of heaters in use decreases as flow rate decreases. The 
necessary voltage variations are obtainable by a General Radio Company 
Type W20HM "Variac" Auto transformer , 0 - 280 volt, 8 amp load from a 240 
volt 50 - 60 cycle line. 

The recommended procedure is to decrease the number of heaters insofar 
as practical while maintaining a high voltage. When there are only four 
heaters in use, voltage reduction becomes necessary to obtain the required 
20° temperature rise, keeping the number of heaters constant. Meeting 
the current requirement is the dictating factor. It is noteworthy that the 
heaters should never be energized unless sufficient air is flowing, other- 
wise the wires will sag and short circuiting will result. 
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MATRIX HOLDER AND TEST SECTION : (See Figures 31, 32, and 35) 



The matrix holder and test section are made of polyethylene plastic. 

The test section casing holds the T^ and T^ thermocouple grids permanently 
mounted, as well as the upstream and downstream static pressure taps required 
for determining friction factor. All parts were machined to close tolerance 
to assure a snug fit for the matrix holder and to guarantee good alignment 
of the flow channel through the heater and test sections. The flow channel 
is 3 - 3/16 M x 3 - 3/16 M and matrices of flow lengths up to 3 M may be 
tested. As mentioned previously under Temperature Measuring System, the T^ 
thermocouple grid is located in the matrix holder. 

Correct positioning of the test matrices in the flow channel is en- 
sured by using styrofoam plastic inserts. In addition to the aforemention- 
ed systems, one item remains to be considered - the Inlet Cone and Screen 
Straightener Section. The inlet cone was designed using Smith and Wang 
[18] as a guide. Although this reference pertained to contracting cones for 
circular duct, the resultant inlet cone for this rig proved very satisfactory. 
To assist in guaranteeing a uniform velocity at the throat, a stainless 
steel 60 mesh screen was installed (See Figure 1). Space is available in 
the screen straightening section to add additional screens if they should 
be required for flow rates beyond the range of this investigation. 
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APPENDIX C 



PERFORATED NICKEL GEOMETRIC PROPERTIES 
The material used in the perforated plate matrices was pure nickel, 
electro-deposited sheet of integral structure manufactured by Perforated 
Products, Incorporated. Types 160/40TV, 160/40Q, 50G, 125M and 125P 
were utilized in the matrices - the first two of these having slotted 
openings and the remaining three having round openings. A brief descrip- 
tion of each of the aforementioned types as specified by the manufacturer 
is set forth in Table C-l. 

Inasmuch as the perforations in all cases were conical, the use of 
the average slot length or average hole diameter was considered a propos. 
Close examination of each of the different types revealed that this 
average value was somewhat less than the dimensions specified by the manu- 
facturer in each case. In addition to the requirements for determining 
porosity of the plate, a correction for the increased longitudinal conduc- 
tion length due to the perforations as well as a correction to the solid 
conduction cross-sectional area was performed. The following idealized 
patterns for each type used, best illustrate these corrections. 

16Q/40TV 
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Cross Hatched Area, A ^ = (.0063) (.0,25) = .0001575 in. 2 

2 

Slotted Area, A gl , in A ch - (.0008) (.0155) - .000012^ in. 
Plate Porosity, the fraction of open or slotted area, 




Solidity, s, 



is 



1 - 



s = 1 




.0000124 

.0001575 



= 0.9212 




, conduction path length, = 0.0125 + 0.0063 = 0.0188" 



L. 



k. 

i 



.0063 

.0188 



= 0.3351 



The solid cross sectional area for heat conduction, A^ varies 

i 

along the path length. The area, A^ , selected was the cross sectional 

i 

area between the slots normal to the flow direction, i.e., 

= (0.0095) (0.0022) = 0.0000209 in. 2 • 
i 

Just as L = 1^ where L is the total matrix flow length and 

is the flow length between two successive perforations; so does 

i 
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and A^ = ^ A^ . In evaluating A^, a conduction area reduction ratio is 

1 , ^ , j> cross sectional area of perforations; 

employed and is defined as ^ = cross sectional area 

therefore, A , is the solid cross sectional area multiplied by 1 

K. « 

i. e. , 



A k - A s (1 - / ) 



= 1.6058 
= 0.00424 ft. 



( 



1 - 
2 



(.0155) (.0022 )1 

(.025) (.0022) j 



= 0.6102 in 



2 



160/40Q 




A . = (.0063) (.025) = .0001575 in. 2 
ch 

2 

A = (.0019) (.0161) * .00003059 in. 
si 

1/ = .1942 
s = .8058 

L = 0.3351, same as 160/40TV 
^ = . 644 

? 

A^ = A (l - ^ ) = .002887 ft. 
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125M 



h- .0079" — #-l 




A = (.00684) (.0079) = .00005404 in. 2 
ch 

A hole = H (.00366) 2 = .00001052 in. 2 

4 

2/ = A hole = .00001052 =0.19466 

A . .00005404 

ch 



s = .8053 



|<- .0079" — J 
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The equivalent conduction length, L. , is computed in the case of 

K i 

round perforations assuming the path to be as shown above, i.e., 2 TT 

““3 



where r^ is the radial distance equal to one-half the center to center 
distance. Thus 



Li, = 2 'TV (.0079 )= .00827" 

k i 3 2 



and T 

L i 

L = .00684 = .82709 
k i .00827 

A is determined as in the case of the slotted perforations, 
k 

= diam. of hole x plate thickness 

center to center distance x plate thickness 



= (.00366) (.0016 ) = .4633 

(.0079) (.0016) 

Therefore, A = A (1 - ^ )= 1.1679in.^ (1 - .4633) 

k s 

= .6268 in . 2 = .00435 ft . 2 



125P 
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A ch = (.00684) (.0079) = .00005404 in 2 

A hole “ ^ (.00291) 2 - .00000665 

4 

Z' = 0.12307 
s = 0.8769 
L ki = .00827" 

L i = .82709 
L ki 

£ = 0.3684 

A = A (1 - ) = .00640 ft 2 

tC s 

50G 



|*_ . 0197 " _^| 
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^ch 

^hole 

TS = 



s = 




(.01706) (.0197) = .00033608 in. 2 

2 2 
IT (.0138) = .00014957 in. 

4 

.44504 

.55496 

.02063" 

.82697 

.7005 

A o (1 - 4 ) = .00243 ft. 2 

S ' 
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Type Opening Hole Diameter Plate Thickness Center to Center 7. Open Area 

_____ (inches ) (inches) Distance, (inches) 
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APPENDIX D 



SAMPLE CALCULATIONS 

The following sample calculations are based on data obtained from 
Run #1 for the modified parallel plate matrix constructed from the 160/40Q 
perforated nickel material. The basic geometric parameters are illus- 
trated in Appendix A. 

Prior to the reduction of data, several fixed parameters must be 
evaluated. 

DIMENSIONS OF FINNED SHEETS: 3.9375" x 2.00" x .0016" 

(BEFORE FORMING) 

DIMENSIONS OF SPLITTER PLATES: 3.1875" x 2.00" x 0016" 

NUMBER OF FINNED SHEETS: 102 

NUMBER OF SPLITTER PLATES: 103 

WEIGHT OF MATRIX: 311.0979gms/453. 6gm/lb = 0.68584 lb. 

MATERIAL CONSTANTS: 

k = 38. 7 Btu/hr - F - ft 
s 

c = 0.1065 Btu/lb - F 
s 

? = 0.321 lb/in 3 = 554.7 lb/ft 3 
s 

PLANE SURFACE AREA, A*, = (102) (2) (2) (3. 9375 ± 3.1875) - 

144 

20.1875 ft 2 

HEAT TRANSFER AREA, A, = PLANE SURFACE AREA x SOLIDITY 

A = A*s 

= (20. 1875) (.8058) 

A = 16.2667 ft 2 

I 2 

FRONTAL AREA, A , = 3.141" x 3.1875" = 10..012 in 
fr i 

2 

A fr = 0.06953 ft 
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It has been shown in Appendix A; however, that this can be alter- 



nately expressed as: 



D H = 4 



Inasmuch as the perforations have no effect on the hydraulic diameter 
the value of @ rather than ^ is used. 



D„ = (4) (0.8834) = 

(1742.1) 



0.002028 ft 



Now, we are in a position to reduce the data. The following presenta- 
tion outlines the steps required; however, as far as practicable the 



actual evaluation is performed by 
in Appendix E. 


the digital computer program set forth 


RECORDED DATA: 


d = 2.310" 


T = 69.0 °F 


o 


o 


AP q - 6.69 in H 2 0 


= 30.140 in Hg. 

atm 


A P = 10.15 in HO 
m 2 


Chart Speed = 4 in/sec. 


P = 10.65 in HO 


(3 = d /d = 2.310 - 0.75 


s 2 


° 3.08 


P Q = 20.80 in H 2 0 


DETERMINATION OF MASS RATE OF FLOW 


: (See Appendix A) 


2 


_ i ' 


rr, - 589.81 Fa Y d 0 C 


1/ AP o 




>| l -$ 4 


If T o 


where: C = C + AC ["" 10 

L %(d o )J 


from reference [l5^ 



P = (P „ - P /13.6) 0.4912 ; T - t (°F) + 459.7 

3 atm o o o 

Y =" 1.0 Fig. 40b of Reference [l6] 

F fl "= 1.0 Fig. 30 of Reference [l6] 



103 



For 

ly 



@ =0.75 and d = 3.08", 



P- = 30.140 - 20.80 

l 13.6 J 



C q => 0.60691 and AC = 0.03839, respective 
0.4912 = 14.0537 psi 



SOLID MATRIX CROSS SECTIONAL AREA, A g , 

A g = 102 x 3.9375" x .0016" + 103 x 3.1875" x .0016" 
= 1.1679 in 2 
A g = 0.00811 ft 2 



FREE FLOW AREA, A , = FRONTAL AREA - SOLID CROSS SECTIONAL AREA 



A = A, - A 
c fr s 



A c = 0.06142 ft' 



CONDUCTION AREA CORRECTED FOR EFFECTS OF PERFORATIONS, A, , 

k 

A k = A (1 - ^ ) where ^ is the conduction area reduction 
s ratio (See Appendix C). 

= 0.00811 (1-0.644) 

A k = 0.00288 ft 2 



POROSITY, p, = MATRIX FLOW VOID VOLUME 

MATRIX VOLUME 



P = 



0.00135 

0.01159 



0.8834 



= Vm - Vs 
Vm 



Vs 

Vm 



COMPACTNESS, (3 



- TOTAL HEAT TRANSFER AREA 
MATRIX VOLUME 



= A 



Vm 



16.2667 

0.01159 



1403.75 ft 2 /ft 3 



Note that (3 includes the effect of area reduction due to perfora- 



tions. 

(3 for an unperforated surface would therefore simply be the plane 
surface area divided by the matrix volume. 
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8 - A* 
v Vm 



20.1875 

0.01159 



1742.1 ft 2 / ft 3 



THE HYDRAULIC DIAMETER, D , = 4r, = 

H h 

4 FLOW CROSS SECTIONAL AREA 
WETTED PERIMETER 



T = 69.0 + 459.7 = 528.7 °R 
o 

m = (589.81 (2.310) 2 C 

1 1 - (0.75) 4 
m = 1605.242 C 




.69 x 14.0537 
528.7 



\ 



assume C = 0.6000 



then "•> = 963.145 lb/hr 
C = C o + A C T 10^ 

L NR < d o). 



where N R ( d ) = (4) (12) m 

(3.1416)(3.08)yU 

= 4.961 tb 

r 






is evaluated at temperature t by means of the linear approximation 




0.0395 + 0.64167 x 10' 4 t = 0.0395 + 0.64167 x 10" 4 (69.0) 

o 



= 0.04392 lb/ft - hr 



therefore, 

C = 0.60691 + 0.03839 



(10 4 ) (0.04392) 
(4. 961) (963. 145) 



C = 0.6104 

This value of C is compared with the assumed value of C and successive 
iterations are performed until the two values are equal. The iterative 
process shows that a value of C = 0.6172 satisfies this requirement and 
thus 

rn = 990.756 lb/hr 
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MATRIX REYNOLDS NUMBERS : 

Two values for Reynolds number are required, one for the isothermal 
flow friction factor evaluation and second for the heat transfer evalua- 
tion. The two values differ only in so far as the difference in the 
absolute viscosities. A is evaluated at orifice temperature (t Q ) 
whereas yU H (heat transfer) is evaluated at average bulk fluid tempera- 
ture which is assumed to be at t + 10 °F. 



n R£ = 4Lm_ 

A /*f 



A Hu 



(4) 12 (990.756) = 

(16. 2667) (.04392) 



924.54 



= N Rf/f!£ 



= 924.54 (.04392 ) = 

(.04456) 



911.23 



CONDUCTION PARAMETER: 






k s A k 
Li^c 



- (38. 7) (0.00288) = 0.00282 
(2/12) (990. 756) (. 24) 



therefore 



X_L_ 



L = 0.3351 from Appendix C. 



Tu- 



(.00282) (.3351) = 



0.000944 



MAXIMUM SLOPE : 

The slope of the cooling curve was determined by the method pres- 
cribed in Appendix A from Figure 36. Slopes were determined for each of 
the curves and the average value was used. 



MAXIMUM SLOPE = C g 

JT 



w 



max 



x SLOPE x CHART 
SPEED 



C = W c 
s s s 



(0.68584 lb) (0.1065 Btu/lb *F) = 0.07304 



C = = (990.756 lb/hr) (0.24 Btu/lb *F) _1 hr_ = 0.06605 



3600 



sec 
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I 




























































SLOPE = 0.14357 (See Figure 36) 



whence, 

MAXIMUM SLOPE = (.07304) (0.14357) (4) 

(.06605) 



0.63507 



COLBURN i FACTOR : 

Enter Table XI with the values of \ and MAXIMUM SLOPE to deter- 

k 

mine N^ = 4.05. Linear interpolation from this table is sufficient. 

j = Nij u A c Np£ 2/3 from Appendix A, 

A 

2/3 
PR 

therefore 



N = 0.8017 - 0.82353 x 10 t 0 using a linear approximation; 



j = (4.05) (0.06142) (.7960) = 

(16.2667) 



0.01217 



FANNING FRICTION FACTOR : 

From Appendix A, the friction factor is: 

f - [***, (J$) - (Kc-Ke) V) 

P at m = 30*140 in Hg x 0.4912 psi = 14.8048 psi 

in Hg 



L, 



P = 10.65 in HO x 5.204 psf = 55.42 psf = 0.3849 psi 
S 2 in H 2 0 



P - 10.15 in H o 0 x 5.204 psf 
m l 



psf 

in H 2 0 



= 52.82 psf = 0.3668 psi 



P, = P „ - P = 14.8048 - 0.3849 = 14.4199 psi 

i atm » 

P 2 = Pj_- AP m = 14.4199 - 0.3668 = 14.0531 psi 

P mean = P 1 + P 2 = 14.2365 psi 

2 

P mean » ( p m ean> < 144 > " M r..?.M5.1Q44) = 0.072681b/ft 3 

- E<t V - 459. 7) (53. 35) (528. 7) 

o 

G = rn = 990.756 = 16130. 841b/hr - ft 2 = 4.481 lb/sec - ft 2 

A 0.06142 

c 
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H - 0.000507 = 0.00304 

r 2/12 



K = 0.48 and K = -0.33 from Figure 38. 

c e 



( 1 + P ) 
f = 



1 + (0.8834) = 1.7804 



(2) (32. 2H0.07268) (52.82) - (0.48 - 0.33) - .3668 

(4.481)* 14.2365 



(1.7804) (0.00304) 

= (12.3126 - 0.150 - 0.04587) (0.00304) 
f = 0.03681] 

The thermal properties, ja . and , were obtained by linear 

interpolation of data from reference j[6^ . The values of k g and c g 
were obtained directly from references j^5^ and ^4~J for nickel and 
stainless steel respectively. 
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APPENDIX E 



FORTRAN 60 DATA REDUCTION COMPUTER PROGRAM FOR A CDC 1604 DIGITAL COMPUTER 

A computer program was formulated in order that as much data reduction 
as practicable would be performed by the digital computer. (See Table E- 
1 ). 

It was hoped that this program together with a subroutine would produce 
Colburn j factor directly. An unsuccessful attempt was made to fill out 
table values of N^ u (TABLE XI) by curve fitting using Aitken's iterative 
method for fitting an nth degree Lagrange polynomial to the n + 1 given 
points - allowing for unequal spacing. As a result of reference [8^J , this 
approach was abandoned; however, the computer program presented is still 
helpful in that values of conduction parameter and maximum slope are end 
results. 

Entering TABLE XI with these parameters yields a value of N^, which 

when multiplied by Cl (Cl ~ Np r ^^ * A^/A) produces the Colburn j factor. 

2/3 

A linear approximation for values of viscosity and N^ as functions 

of temperature is included. An added convenience in the computer program 
is the evaluation of the orifice constant, K. Using ^15^ , an assumed 
value for C (a slowly varying function) was made. The correct value for 
C was then evaluated by an iterative process by the computer. With the 
value of C determined, K (a rapidly varying function) was found from the 
following simple relationship: - _ C 

where Q is the ratio of the orifice diameter to the inside pipe dia- 
meter. 

With the computer program presented, an unlimited number of matrices 
can be reduced by simply changing the first card of the program. (L = 
number of matrices for which data is to be reduced). 
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The following is a list of the input data for the FORTRAN 60 program: 
(See sample input data in TABLE E - 2) 

CARD NO. 



30 XML 


Matrix length, ft 


AS 


2 

Solid conduction area, ft 


RH 


Hydraulic radius, ft 


POR 


Porosity 


AFL 


2 

Free flow area, ft 


AHT 


2 

Matrix heat transfer area, ft 


WM 


Matrix weight, lbs. 


40 EK 


Exit flow coefficient 


CK 


Entrance flow coefficient 


SK 


Solid thermal conductivity, Btu/hr ft 0 F 


CM 


Matrix material specific heat, Btu/lb °F 


XLR 


Ratio of L/L^ (Flow length/Equivalent conduc- 
tion length) 


N 


Number of runs 


10 DO 

• 


Diameter of orifice plate opening, inches 


DELPO 


Pressure drop across orifice, inches H^O 


P0 


Static pressure upstream of orifice, inches ^0 


ATMP 


Atmospheric pressure, inches Hg 


BETA 


Ratio of orifice diameter to inside pipe diameter 


10+2 CS 


Chart speed, sec/inch 


DELH 


Pressure drop across matrix, inches ^0 


HS 


Static pressure at matrix inlet, inches H^0 


SLO 


Slope of generalized cooling curve, (inches)"* 


TEMPO 


Fluid temperature at orifice (t Q ), °F 
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The output from the program is as follows: 



EMDOT 


Mass rate of flow (m), lb/hr 


RNUMP 


Reynolds number for the pipe 


RNUMH 


Reynolds number for heat transfer 


CAPFL 


Thermal capacitance of the fluid (C), Btu/sec °F 


CAPS 


Thermal capacitance of the solid ((T s ), Btu/ °F 


CONDPAR 


Conduction parameter ( A. ) 


SLOMAX 


Maximum slope of generalized cooling curve 


CONPARE 


Equivalent conduction parameter ( A. K ) 


G 


Mass flow velocity (m/A c ), lb/sec - ft L 


RNUMFR 


Reynolds number for friction 


FFR 


Fanning friction factor 


ORIFICE K 


Orifice constant, K 


PRANDTL NO. 


2/3 (Prandtl number) 2/ , 3 N / 3 


Cl 


2/3 

Parameter s N x A„/A 

Pr c 


UFR 


Fluid viscosity, f*- , lb/hr-ft 

(evaluated at fluid temperature at orifice) 


ESTD 


Flow friction power per unit area (E __), HP/ft 2 

Oil/ 


STDH 


Parameter for evaluating heat transfer power per unit 
area (h STD ), Btu/hr ft 2 °F. (STDH x N Tu = h STD ) 
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PROGRAM SINGBLO 
L = 1 . 

K= 1 

30 READ 300*XML*AS»RH» POR * AFL » AHT > WM 
300 FORMAT (7F 10. 5) 

40 READ 40 0 » E K » C K » S K » CM » X L R » N 
400 FORMAT ( 5F 10. 5, 12 ) 

M= 1 

10 READ 10G,DO,DELPO,PO,ATMP»3ETA 

100 FORMAT ( 5F10 .3 ) 

READ 101»CS>DELH , HS , SLO , TEMPO 

101 FORMAT ( 3F 10 . 5 ) 

P3= ( ATMP- ( PO/ 13 .6 ) ) * . 49 1 2 
TU=TEMPO+459»7 
UFR = 0.0395 + 0 . 641 67 E-4*TEMP0 
TB = T EMPO+ 10 • 0 
UH=0. 0395+0. 64167E-4+TB 
WRITE OUTPUT TAPE 4,9 
9 FORMAT ( 1H1 ) 

WRITE OUTPUT 'TAPE 4,4,M 
4 FORMAT (9H RUN NO. 12///) 

11 IF(AB5F (BETA-. 75)-. 000001)2 7, 27, 12 

12 IF (ABSF (BETA-. 64)-. 000001 ) 26,26, 13 s 

13 IF (ABSF (BETA-. 50) -.00 0001 ) 2 5-, 2 5 , 14 

14 IF(ABSF(bETA-.4C)-.0u0001 124,24,15 

15 IF (ABSF (BETA-. 25) -.000001) 23, 2 3, 16 

16 IF (ABSF (BETA-. 15) -.000001 >22,22,17 
17. GO TO 8 0 

21 X2=0. 59225 
DELX = O'. 01916 
ASSUM=0 .6000 
GO TO 2.8 

22 X2=0. 59171 
DELX=0. 01691 
ASSUM=0 .6000 
GO TO 28 

23 X2=0. 59184 
DELX=0. 01622 
ASSUM= 0.6000 
GO TO 28 

24 X2=0. 59448 
DELX=0. 01867 
ASSUM= 0.6000 

• GO TO 28 



TABLE E-l 



112 



25 X2=0. 59850 
DELX=0. 02125 
ASSUM* 0.6000 
GO TO 2e 

26 X2=0. 60575 
DELX=0. 02595 
ASSUM=0.60C0 
GO TO 28 

27 X2=0. 60691 
DELX=0. 03839 
ASSUM= 0.6000 

28 COR=ASSUM 
NC = 0 

B1 =1 ../SQRTF ( l.-BETA**4> 

31 NC=NC+ 1 

IF (NC-4)29»29>60 

29 EMDOT=589.81*DO*DO*SORTF( DELPO*P3/T0 )*COR*bl 
RNUMO= ( 4.961 *EMDOT/UFR ) /SETA 
CORA=X2+DELX*SQRTF( 10000. /RNUMO ) 

WRITE OUTPUT TAPE 4, 8, CORA 
8 FORMAT (6H COR A= F 1 0 . 5 / / ) 

32 COR=CORA 
GO TO 31 

60 PR=0.8017-. 82353E-4*TEKPO 
C1=PR*AFL/AHT 
CP = 0 . 2 4 

RNUMP = 5 . 09* EMDOT/UFR 

RNUMH= ( 4 . 0*XML* EMDOT ) / ( AHT *UH ) 

RNUMFR=RNUMH*UH/UFR 

CAPFL=EMDOT*CP/ 3600.0 

CAPS=WM*CM 

SLOMAX= ( CAPS*SLO) / (CAPFL*CS ) 

CONDPAR = ( SK*AS ) / ( XML*EMDOT*CP ) 

CONPARE=CONDPAR*XLR 

3=CK+EK 

C=1 .0+POR**2 

G=EMDOT / ( AFL*3600.0 ) 

DELP = DELH*.036 13 
PS=HS*. 03613 • 

PA=ATMP* . 49 1 2 
P1=PA-PS 
P2=P1-DELP 
PV= ( P1+P2 ) /2. 

RHOM=( PM* 144.0) / ( 53. 3* ( TEMPO + 460.0 ) ) 



TABLE E-l 
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FF 1 = ( (64.4*RH0M*DELP* 144.0 )/(G*G) ) 

FF2=C*DELP/PV 
FF3=FF1-FF2-B 
FFR=FF3*RH/XML 
CONST = COR*B 1 

STDH=(0.02195*RNUMH*C1 >/(4.*RH) 

ESTD= 1 ■ HE-7*( 1.0/(4.*RH) ) **3*FFR* ( RNUMFR/ 1000 . )**3 
WRITE OUTPUT TAPE 4, 7 , STDH , ESTD 
7 FORMAT ( 1 9H CONSTANT FOR HSTD=F 10 . 5 » 3X » SH ESTD=F10 . 5/// ) 

WRITE OUTPUT TAPE 4 , 1 , EMDOT , RNUMP , RNUMH , CAPF L , CAPS 

1 FORMAT ( 7H EMD0T=F10. 5 » 3X,7H RNUMP=E20 . 8 , 3X , 7H RNUMH=F 10 . 5 , 3X , 

17H CAPFL=F1C.5»3X»6H CAPS=F10. 5/// ) 

WRITE OUTPUT TAPE 4 , 2 * CONDPAR , SLOMAX , COMPARE 

2 * FORMAT ( 9H C0NDPAR=F10.5»3X»8H SLOMAX = F 10 . 5 , 3X , 9H CONPARE = F 1 0 . 5/ / / ) 

WRITE OUTPUT TAPE 4 , 3 , G , RNUMFR , FFR , CONST 
30 FORMAT ( 3H G = F 10 . 5 » 3X » 8ri RNUMFR = F 10 . 5 , 3X , 5H FFR=F10.5»3X> 11H ORIFIC 
IE K=F10.5///) 

WRITE OUTPUT TAPE 4 , 5 , PR , C 1 , UFR 

5 FORMAT ( 16H PRANDTL NO . 2 / 3 = F 1 0 . 5 , 3X , 4H C 1 =F 10 . 3 » 3X , 5H UFR = F10 . 5/// > 
M=M+1 

IF(N-M) 50,10,10 
50 K=X+ 1 

I F ( L-K ).20 , 30 , 30 

80 PRINT 81 

81 FORMAT (30H ORIFICE CONSTANT ERROR EXISTS) 

GO TO 20 

20 STOP 
END 

END < 



TABLE E-l 



114 



SAMPLE INPUT DATA FOR PROGRAM SINGBLO 



rs 



0.1667 


.00640 


.00049 


.8542 


.05939 


17.70242 


-.32 


.51 


38.7 


.1065 


.82709 


11 ' 


2.310 


6.78 


21.5 


30.140 


.75 




.25 


9.30 


12.2 


.17768 


68.5 




2.310 


4.40 


15.45 


30.120 


.75 


' 


.25 


7.15 


8.30 


.15259 


71.0 


' 


2.310 


2.74 


10.95 


30.120 


.75 




.25 


5.15 


5.80 


.12882 


73.0 




1.971 


3.82 


7.10 


30.120 


.64 




.50 


3.40 


3.50 


.21040 


72.0 




1.540 


6.89 


4.50 


30.120 


.50 




. 50 


2.33 


2.06 


.17379 


72.7 




1.540 


2.91 


2.60 


30.110 


.50 




.50 


1.42 


0.95 


.14039 


72.7 




1.232 


2.75 


1.40 


30.110 


'.40 




.50 


0.80 


0.43 


.09232 


72.5 




0.775 


6.00 • 


0.68 


30.110 


.25 




1.00 


0.41 


0.18 


.10170 


69.0 




0.775 


5.98 


0.67 


30.110 


.25 




1.00 


0.40 


0. 178 


.10316 


69.0 




0.775 


1.78 


0.35 


30.110 


.25 




2.00 


0.22 


0.075 


.10321 


70.0 




0.462 


2.60 


0.15 


30.110 


.15 




4.00 


0.080 


0.025 


• 07513 


70.0 






t 

\ 


TABLE 


E - 2 




1 







.4844 
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